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ABSTRACT 

Whether gamma-ray bursts are highly beamed or not is a very difficult but 
important problem that we are confronted with. Some theorists suggest that 
beaming effect usually leads to a sharp break in the afterglow light curve during 
the ultra-relativistic phase, with the breaking point determined by 7 = l/#o, 
where 7 is the Lorentz factor of the blastwave and 8q is the initial half opening 
angle of the ejecta, but numerical studies tend to reject the suggestion. We 
note that previous studies are uniformly based on dynamics that is not proper 
for non-relativistic blastwaves. Here we investigate the problem in more detail, 
paying special attention to the transition from the ultra-relativistic phase to 
the non-relativistic phase. Due to some crucial refinements in the dynamics, 
we can follow the overall evolution of a realistic jet till its velocity is as small 
as f3c ~ 10~ 3 c. We find no obvious break in the optical light curve during the 
relativistic phase itself. However, an obvious break does appear at the transition 
from the relativistic phase to the Newtonian phase if the physical parameters 
involved are properly assumed. Generally speaking, the Newtonian phase is 
characterized by a sharp decay of optical afterglows, with the power law timing 
index a ~ 1.8 -- 2.1. This is due to the quick lateral expansion at this stage. 
The quick decay of optical afterglows from GRB 970228, 980326, and 980519, 
and the breaks in the optical light curves of GRB 990123 and 990510 may 
indicate the presence of highly collimated 7-ray burst ejecta. 

Subject headings: gamma rays: bursts - - ISM: jets and outflows - 
hydrodynamics — radiation mechanisms: non-thermal — relativity 



-3- 

1. Introduction 

Till the end of August 1999, X-ray, optical, and radio afterglows have been observed 
from about 16, 11, and 5 gamma-ray bursts (GRBs) respectively (Costa et al. 1997; Bloom 
et al. 1998; Groot et al. 1998; Kulkarni et al. 1998, 1999; Harrison et al. 1999; Stanek et al. 
1999; Fruchter et al. 1999; Galama et al. 1999a). The so called fireball model (Goodman 
1986; Paczyhski 1986; Meszaros & Rees 1992; Rees & Meszaros 1992, 1994; Katz 1994; 
Sari, Narayan, & Piran 1996; Dermer & Mitman 1999; Dermer, Bottcher, & Chiang 1999, 
2000; Dermer 2000) is strongly favored, which is found successful at explaining the major 
features of GRB afterglows (Meszaros & Rees 1997; Vietri 1997; Tavani 1997; Waxman 
1997a; Wijers, Rees, & Meszaros 1997; Sari 1997a; Huang et al. 1998a, b, 1999a, b; Dai 
& Lu 1998a, b, c; Dermer, Chiang, & Bottcher 1999; Dermer, Chiang, & Mitman 2000). 
However, we are still far from resolving the puzzle of GRBs, because their "inner engines" 
are well hidden from direct afterglow observations. 

To unveil the nature of the "inner engine", we first need to know the energetics 
involved in a typical burst, which itself depends on two factors: (i) the distance scale 
of GRBs, this has been settled since the BeppoSAX discovery of GRB 970228; (ii) the 
opening angle of GRB ejecta, i.e., whether gamma-rays are radiated isotropically or not, 
this question is still largely uncertain. Most GRBs localized by BeppoSAX have indicated 
isotropic energies of 10 51 -- 10 52 ergs, well within the energy output from compact stellar 
objects of solar-mass. However, GRB 971214, 980703, 990123, and 990510 have implied 
isotropic gamma-ray releases of 3.0 x 10 53 ergs (0.17 M & c 2 , Kulkarni et al. 1998), 1.0 x 10 53 
ergs (0.06 M Q c 2 , Bloom et al. 1998), 3.4 x 10 54 ergs (1.9 M G c 2 , Kulkarni et al. 1999; 
Andersen et al. 1999), and 2.9 x 10 53 ergs (0.16 M Q c 2 , Harrison et al. 1999) respectively. 
Moreover, if really located at a redshift of z > 5 as suggested by Reichart et al. (1999), 
GRB 980329 would imply an isotropic gamma-ray energy of 5 x 10 54 ergs (2.79 M Q c 2 ). 
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Such enormous energetics has forced some theorists to deduce that GRB radiation must 
be highly collimated in these cases, with half opening angle 9 < 0.2, so that the intrinsic 
gamma-ray energy could be reduced by a factor of 10 2 -- 10 3 , and could still come from 
compact stellar objects (Pugliese, Falcke, & Biermann 1999). Obviously, whether GRBs are 
beamed or not is of fundamental importance to our understanding of their nature. 

How can we tell a jet from an isotropic fireball? Gruzinov (1999) has argued that 
optical afterglows from a jet can be strongly polarized, in principle up to tens of percents, 
if co-moving magnetic fields parallel and perpendicular to the jet have different strengths 
and if we observe at the right time from the right viewing angle. More direct clues may 
come from afterglow light curves. Rhoads (1997, 1999a, b) has shown that the lateral 
expansion (at sound speed) of a relativistic jet (with a Lorentz factor 7 > 2) will cause the 
blastwave to decelerate more quickly, leading to a sharp break in the afterglow light curve. 
The power law decay indices of afterglows from GRB 980326 and 980519 are anomalously 
large, a ~ 2.0 (Groot et al. 1998; Owens et al. 1998), and optical light curves of GRB 
990123 and 990510 even show obvious steepening at observing time t > 1 -- 2 d (Kulkarni 
et al. 1999; Harrison et al. 1999; Castro-Tirado et al. 1999). Recently GRB 970228 was 
also reported to have a large index of a ~ 1.73 (Galama et al. 1999b). These phenomena 
have been widely regarded as the evidence of the presence of relativistic jets (Sari, Piran, & 
Halpern 1999; Castro-Tirado et al. 1999). 

However, numerical studies of some other authors (Panaitescu & Meszaros 1998; 
Moderski, Sikora, & Bulik 2000) have shown that due to the increased swept-up matter and 
the time delay of the large angle emission, the sideway expansion of the jet does not lead 
to an obvious dimming of the afterglow. Thus there are two opposite conclusions about the 
jet effect: the analytical solution predicts a sharp break, while the numerical calculation 
shows no such sharp breaks. It is very interesting to note that a recent analytic treatment 
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by Wei & Lu (2000) shows that the sharp break predicted by Rhoads is usually not exist 
unless the beaming angle is very small. This analysis seems to have given some supports to 
the numerical studies. 

We find that previous studies on jet effects need to be improved urgently in the 
following three aspects: (i) Afterglows during the mildly relativistic (2 < 7 < 5) and 
non-relativistic (7 < 1.5) phases of the expansion are of great importance to us, since they 
may correspond to observing time of£~2--5d (Huang et al. 1998a, b). However, 
conventional dynamical model can not transit correctly from the ultra-relativistic phase 
(7 ^> 1) to the non-relativistic phase for adiabatic shocks. This has been stressed by Huang, 
Dai, & Lu (1999, 2000). Therefore previous numerical studies were based on improper 
dynamical equations. They could describe an ultra-relativistic jet, but they gave spurious 
results in the mildly relativistic and non-relativistic phases (Huang, Dai, & Lu 1999a, b); 
(ii) It is reasonable to assume that the lateral expansion speed of the beamed ejecta is 
just the co-moving sound speed c s . Usually we take c s = c/v3, where c is the speed of 
light (Rhoads 1997, 1999a). However in realistic case we expect c s to vary with time, and 
especially it will by no means be c/\^3 when 7 ~ 1. This is another reason that previous 
studies are not proper for mildly relativistic and non-relativistic jets; (iii) In previous 
studies, the expansion of the beamed ejecta is supposed to be solely adiabatic. In fact the 
blastwave should be highly radiative at first. It evolves slowly into an adiabatic blastwave 
at later stages (Dai, Huang, & Lu 1999). 

This paper is aiming to overcome the problems mentioned just above. We present our 
refined dynamical model in Section 2. It can describe the overall evolution of a realistic, 
jetted GRB ejecta. Synchrotron radiation from shock accelerated electrons is formulated in 
Section 3. Our detailed numerical results are presented in Section 4. Section 5 is a brief 
discussion. 
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2. Dynamics 

2.1. Basic Equations 

Let R be the radial coordinate in the burster frame; %,t co , and t be the time from 
the event measured in the burster frame, co-moving ejecta frame, and terrestrial observer's 
frame respectively; 70 and M e j be the initial Lorentz factor and ejecta mass and 9 the half 
opening angle of the ejecta. The burst energy is Eq = 7oM e jC 2 . 

It is useful to write down the following simple relations at first: 

dR = (3cdt h , (1) 



dt h = -fdt co = 7(7 + yV - l)dt, (2) 

where (3 = a/7 2 — 1/7. 

The evolution of radius (R) and swept-up mass (to) is described by 

d R 

_ = / g C7 (7 + v VTI )? (3) 

^- = 2irR 2 (l- cos d)nm p , (4) 

dR 

where n is number density of surrounding interstellar medium (ISM) and m p is mass of 
proton. The opening angle 6 is usually evaluated as 9 rs a/R, where a is co-moving lateral 
radius of the ejecta (Rhoads 1999a; Moderski, Sikora, & Bulik 2000). This is not ideal in 
numerical evaluations since 9 may be as large as 0.5 — 0.8 at later times. Here we use the 
following differential equation: 

d9 _ 1 da c s (7 + V7 2 - 1) 



dt Rdt R 



(5) 



As for the evolution of the Lorentz factor 7, the following equation has been widely 

used for both isotropic fireballs and beamed ejecta: 

d'j 7 2 — 1 
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where M is the total mass in the co-moving frame, including internal energy U, 
dM = [(1 - e)7 + e]dm (Chiang & Dermer 1999; Piran 1999; Moderski, Sikora, & Bulik et 
al. 2000). Here e is the radiative efficiency defined as the fraction of the shock-generated 
thermal energy (in the co-moving frame) that is radiated, e = 1 corresponds to highly 
radiative case and e = corresponds to adiabatic expansion. However Huang, Dai, & Lu 
(1999a, b) have demonstrated clearly that for adiabatic expansion equation (^) does not 
agree with the Sedov solution in the non-relativistic limit: equation (^) gives /3 oc R~ 3 while 
the Sedov solution requires (3 oc R~ 3 ^ 2 (Sedov 1969). They have proposed a new equation. 
Now we briefly repeat their derivation. 

In the observer's frame, since the total kinetic energy of the fireball is 
E\i = (7 — l)(M e j + m)c 2 + (1 — t)^U (Panaitescu, Meszaros, & Rees 1998), and 
the radiated thermal energy is €7(7 — l)dmc 2 (Blandford & McKee 1976), we have 

d[(7 - l)(M ej + m)c 2 + (1 - e)-fU] = -67(7 - l)dmc 2 . (7) 

For the item U, it is usually assumed that dU = (7 — l)dmc 2 (Panaitescu, Meszaros, & 
Rees 1998). Equation (|]) has been derived just in this way. Huang, Dai, & Lu (1999a, b) 
suggested that U — (7 — l)mc 2 is a better approximation that makes equation (0) suitable 
for both ultra-relativistic and Newtonian shocks. Using this expression, it is easy to obtain 
(Huang, Dai, & Lu 1999a, b) 

dm M e j + em + 2(1 — e)7m 

Huang, Dai, & Lu (1999a, b) have shown that this equation is acceptable for both radiative 
and adiabatic fireballs, and in both ultra-relativistic and non-relativistic phases. Here we 
use this equation to follow the overall evolution of beamed GRB ejecta. 

Equations ([]), (||), (0), and M) present a thorough description of jet evolution. But 
before evaluating them numerically, we should give the expression for c s and e. 



2.2. Sound Speed 

The lateral expansion is determined by the co-moving sound speed. The simple 
assumption of c s = c/v3 is not proper for our purpose in this paper. We must derive c s 
from c 2 = (dp'/de') s , where p' and e' are co-moving pressure and energy density respectively. 
Kirk & Duffy (1999) have derived 



C s 



2 _ 1V_ 
P' 



(7 - iy 



(9) 



(7- l)p' + 7p'_ 

where p' is co-moving mass density, and 7 is the adiabatic index. Dai, Huang, & Lu 
(1999) gave a simple and useful expression for 7, 7 = (47 + l)/(37). Since e' = 7p'c 2 and 
p' = (7 — l)(e' — p'c 2 ), it is easy to get 

c; = w -i)( 7 -i)__l__ c s. no) 



In the ultra-relativistic limit (7 > 1,7 « 4/3), equation ( |10D gives c 2 = c 2 /3; and in 



the non-relativistic limit (7 ~ 1,7 w 5/3), we simply get c 2 = 5/5 2 c 2 /9. So equation fllOD is 
a reasonable expression and will be used in our model. 

2.3. Radiative Efficiency 

As usual we assume that the magnetic energy density in the co-moving frame is a 
fraction £| of the total thermal energy density (Dai, Huang, & Lu 1999) 

f^^(7-lWc 2 , (11) 

and that the shock accelerated electrons carry a fraction £ c of the proton energy. Here £ e 
is a parameter characterizing the efficiency of energy transport from protons to electrons. 
In realistic case, it is probable that £ e may vary with time, but little is known about the 
detailed mechanism. So we will take £ e as a constant throughout this article. This implies 
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that the minimum Lorentz factor of the random motion of electrons in the co-moving frame 

is 

. mr,(p — 2) , . 

7e,min = fc 7 " 1 ( , + 1, 12 

m c (p- 1) 
where m e is electron mass and p is the index characterizing the power law energy 
distribution of electrons. We here consider only synchrotron emission from these electrons, 
and neglect the contribution of inverse Compton emission because the latter is unimportant 
particularly at late times (Waxman 1997a). The energy of a typical electron is lost due 
to both synchrotron radiation and expansion of the ejecta, thus the radiative efficiency of 
this single electron is given by t'~^/{t'~^ + 1'^ 1 ) (Dai, Huang, & Lu 1999), where t' syn is the 
synchrotron cooling time, t' syQ = 67rm c c/(crT-B /2 7 e ,mm), with o"t the Thompson cross section, 
and t' ex is the co-moving frame expansion time, t' ex = R/^c). Since all the shock accelerated 
electrons carry only a fraction £ e of the internal energy, the radiative efficiency of the total 
ejecta can be given by (Dai, Huang, & Lu 1999) 

e = ^rW (13) 

syn ' ^ex 

For the highly radiative expansion, ^ e » 1 and tg yn <^C t' cx , we have e ~ 1. The early 
evolution of the ejecta is likely in this regime. For an adiabatic expansion, ^ e < 1 or 
t' syn ^> t' cx , we get e ~ 0. The late evolution is believed to be in this regime. So one expects 
that in realistic case the radiative efficiency evolves from about 1 to 0. In this paper we call 
the jet whose radiative efficiency evolves according to equation (|T3| ) a "realistic" one (Dai, 
Huang, & Lu 1999). 



3. Synchrotron Radiation 
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3.1. Electron Distribution 

In the absence of radiation loss, the distribution of the shock accelerated electrons 
behind the blastwave is usually assumed to be a power law function of electron energy, 

dN' 

— - OC 7 C P , (7e,min < 7c < 7c,max), (14) 

d7 e 
where 7 c ,max is the maximum Lorentz factor, 7 e)max = 10 8 (.E>'/1G) _1 / 2 (Dai, Huang, & Lu 
1999), and p usually varies between 2 and 3. However, radiation loss may play an important 
role in the process. Electrons with different Lorentz factors have different radiation 
efficiencies. Sari, Piran, & Narayan (1998) have derived an equation for the critical electron 
Lorentz factor, 7 C , above which synchrotron radiation is significant, 

6nm c c 

7c = ^mt (15) 

Electrons with Lorentz factors below 7 C are adiabatic ones, and electrons above 7 C are 
highly radiative. 

In the presence of steady injection of electrons accelerated by the shock, the distribution 
of radiative electrons becomes another power law function with an index of p + 1 (Rybicki & 
Lightman 1979), but the distribution of adiabatic electrons is unchanged. Then the actual 
distribution should be given according to the following cases (Dai, Huang, & Lu 1999): 

(i) For 7 C < 7 Cjmin , 

(7,- / -(P+ 1 ) (ry . <iy <iy \ (IQ) 

w l /e 5 V /e,min _ ;e — /c,max/ j \ ±KJ J 



die 



Ci= p V _ _ p AT ele , (17) 

/c,min |c,max 



where iV e i e is the total number of radiating electrons involved. 

(fi) For 7e jmin < 7 C < 7c,max, 

dK_\ C ^c P , (7c,min < 7c < 7c) , 

^ lC 3 7- (P+1) , (7c<7c<7c, m ax), 
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where 



C 2 = C 3 / 7 c 



C, 



/c,rmn /c <c /c,max 



(iii) If 7 C > 7c, max , then 



where 



, -l 



7c(p - 1) 



P 



N, 



clc 



dK 



L'47e J wc,min _ 7c _ 7c, max/ j 



^4 - -p^ ir^-^ele- 

7c,min 7e,max 



(19) 
(20) 

(21) 
(22) 



3.2. Relativistic Transformations 

In the co-moving frame, synchrotron radiation power at frequency v' from electrons is 
given by (Rybicki & Lightman 1979) 

.,. ,, V3e 3 B' [V- (dN' c 



P'W 



m c c H 



v 



die 



F - \d le 



(23) 



where e is electron charge, v' c = S^eB' j (4irm c c) , and 



r+oo 

F(x) — x j K 5 / 3 (k)dk, 



(24) 



with K 5 / 3 (k) being the Bessel function. We assume that this power is radiated isotropically, 



dP'(u') P'(is') 



d£l' 



4tt 



(25) 



Let be the angle between the velocity of emitting material and the line of sight and 
define \x = cos©, we can derive the angular distribution of power in the observer's frame 
(Rybicki & Lightman 1979), 

dP(u) 1 dP\v') 1 P'{v') 



dtt 



7 3 (l- / 3/i) 3 dQ' 7 3 (l-/3/i) 3 4vr 



(26) 
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" = W^WY (27) 



Then the observed flux density at frequency v is 

1 (dP(u) A\ 1 1 . . . 

where A is the area of our detector and D^ is the luminosity distance. 



3.3. Equal Arrival Time Surfaces 

Photons received by the detector at a particular time t are not emitted simultaneously 
in the burster frame. This effect has been emphasized by a number of authors (Waxman 
1997b; Sari 1997b; Panaitescu & Meszaros 1998) and may be of great importance to jet 
radiation (Moderski, Sikora, & Bulik 2000; Panaitescu & Meszaros 1999). In order to 
calculate observed flux densities, we should integrate over the equal arrival time surface 
determined by 

t = f ~ — dR = const, (29) 

within the jet boundaries. 



4. Numerical Results 

For convenience, let us define the following initial values or parameters as a set of 
"standard" parameters: initial energy per solid angle E /Q = 10 54 ergs/4vT, 70 = 300 (i.e., 
initial ejecta mass per solid angle M e j/fi = 0.001867M Q /47r), n — 1 cm -3 , £g = 0.01, 
p = 2.5, -Dl = 1-0 x 10 6 kpc, C, c = 1, 6q = 0.2. For simplicity, we first assume that 
the expansion is completely adiabatic all the time (i.e. e = 0, we call it an "ideal" jet, 
distinguishing it from the "realistic" jet defined in Section 2.3). 
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Figure 1 shows the evolution of the Lorentz factor. We see that the ultra-relativistic 
phase lasts only for ~ 10 5 s, this is the period during which equation (|6]) can be safely 
applied. The mildly relativistic phase lasts from ~ 10 5 s to ~ 10 6 ' 5 s. In short, the ejecta 
will cease to be highly relativistic at time t ~ 10 5 - - 10 6 s. This clearly demonstrates the 
necessity to replace equation @ with our improved expression of equation (||). Figure 1 also 
indicates that the transition from the relativistic phase to the non-relativistic phase occurs 
roughly at t ~ 10 6 ' 5 -- 10 7 s. Figure 2 illustrates the time dependence of shock radius. 

EDITOR: PLACE FIGURE Q HERE. 

EDITOR: PLACE FIGURE | HERE. 

In Figure 3 we present the evolution of the jet opening angle 9. During the 
ultra-relativistic phase (t < 10 5 s), 9 increases only slightly. But at the Newtonian stage 
(t > 10 7 s), the increase of 9 is very quick. Again this figure indicates that the critical point 
is roughly at t ~ 10 6 ' 5 s. The usual approximation of 9 ~ a/R is not a good approach 
and our equation (H) is obviously more reasonable. Figure 4 illustrates jet evolution on 
the y - z plane schematically, where z- axis is just the symmetry axis of the jet and the 
lateral expansion is approximately at y direction. Note that the jet is non-relativistic when 
z ~ 3.2 x 10 18 cm. Again we see the quick lateral expansion at the Newtonian stage. In 
Figure 5 we show some examples of equal arrival time surfaces on the y - z plane. Here we 
assume that the angle between jet symmetry axis and the line of sight is # b s = 0. 

EDITOR: PLACE FIGURE | HERE. 
EDITOR: PLACE FIGURE H HERE. 
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EDITOR: PLACE FIGURE | HERE. 

Figures 6 and 7 illustrate the time dependence of radiative efficiency (e) and total 
kinetic energy (E^) for "realistic" jets respectively. Although the highly radiative phase 
usually lasts for very short period (t < 10 4 -- 10 5 s), it dissipates the total kinetic energy 
substantially (up to 90% in less than 10 3 s). This effect may account for the confusing 
phenomena observed in some GRBs: 7-ray energy released in the initial GRB phase (E 7 ) 
is comparable to or even higher than the total kinetic energy inferred from afterglow 
observations (Wijers & Galama 1999). 

EDITOR: PLACE FIGURE | HERE. 

EDITOR: PLACE FIGURE HERE. 

Figure 8 shows the optical light curves, computed for R band observation. In this 
figure, the thick solid line corresponds to an "ideal" jet with "standard" parameters, and 
viewing angle # bs = 0. For comparison, other lines are drawn with only one parameter 
altered or only one condition changed. Please note that at end point of each curve, the 
average electron Lorentz factor is already as small as 7 e ~ 5, corresponding to bulk Lorentz 
factor of 7 ~ 1.002, thus completely in the Newtonian regime. 

EDITOR: PLACE FIGURE | HERE. 

^From Figure 8, we see that in no case could we observe the theoretically predicted light 
curve steepening during the relativistic stage itself (i.e., when t < 10 6 s), consistent with 
previous numerical studies (Panaitescu & Meszaros 1999; Moderski, Sikora, & Bulik 2000). 
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This has been attributed to the effects of equal arrival time surfaces and more and more 
swept-up material. Here we would like to propose another more reasonable explanation. 
Rhoads is correct in saying that the lateral expansion begins to take effect when 7 ~ l/#o- 
In our calculations this occurs at t ~ 10 5 ' 5 s. However, the blastwave is already in its mildly 
relativistic phase at that moment and it will become non-relativistic soon after that (i.e., 
when t > 10 65 s, see Figure 1). So it is not surprising that we could not see any obvious 
breaks during the relativistic phase, they just do not have time to emerge. 

Fortunately we have pointed out that the lateral expansion is even more notable 
in the non-relativistic phase, so we expect that a break should occur in the light curve 
between the relativistic stage and the non-relativistic stage. The dash-dotted line in Figure 
8 proves our deduction. This line corresponds to an "ideal" jet with £ e = 0.1 and with other 
parameters equal to those in the "standard" set. Here the decay in the relativistic phase 
is Sr oc t~ 161 , in good agreement with theoretical results for isotropic fireballs (Wijers, 
Rees, & Meszaros 1997; Sari et al. 1998), and the decay in the non-relativistic phase is 
Sr oc t -214 , also consistent with Rhoads' prediction for lateral expansion effect (but made 
only for relativistic jets). Previous numerical studies are based on dynamical equations 
improper for non-relativistic blastwaves. They could not correctly reveal the break between 
relativistic and Newtonian stages (Panaitescu & Meszaras 1999; Moderski, Sikora, & Bulik 
2000). 

One may ask why other lines in Figure 8 do not show any breaks. The reason is they 
all correspond to £ e — 1 an d their peaks appear at relatively late stages. For example, the 
thick solid line peaks at about 10 5 s. Then we can not see the initial power law decay (i.e., 
a ~ 1.1) in the relativistic phase and can only see the fast decay of Sr oc t -214 at lage 
stages. However, we should note that although they do not show any obvious breaks, they 
all have a common characteristic: a large timing index of a ~ 1.8 -- 2.1. We suggest that 
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the sharp decline of afterglows from some GRBs (with a ~ 1.7 — 2.0) itself may just be the 
evidence of the presence of a jet (Huang, Dai, & Lu 2000a, b, c). 

Dai & Lu (1999, 2000) have discussed the optical afterglows from isotropic fireballs 
during the non-relativistic phase extensively. The most obvious difference between their 
results and ours is that our timing index at non-relativistic stage is reasonably larger. 

5. Conclusion and Discussion 

In this paper we investigate the detailed dynamical evolution of jets and their 
afterglows. Our model is simple in form and is easy for numerical evaluations. Comparing 
with previous studies, the model is refined in the following aspects: 

(i) Equation (|6]) has been widely used in previous studies, however, it is not correct for 
non-relativistic ejecta (Huang, Dai, & Lu 1999a, b). The dynamics here (i.e., mainly 
our equation (|8|) ) is applicable to both ultra-relativistic and Newtonian jets, so we 
could follow the overall evolution (from 7 ^> 1 to 7 ~ 1) of a jet by using a single set 
of differential equations. Numerical results indicate that the ejecta will cease to be 
ultra-relativistic 10 5 -- 10 6 s after the main GRB. We should consider this refinement 
seriously. 

(ii) We describe the lateral expansion of jets with a refined sound speed expression, which 
gives reasonable approximations during both ultra-relativistic and Newtonian phases 
(Kirk & Duffy 1999). 

(iii) The remnant here is more "realistic", i.e., the radiative efficiency evolves from 1 
(corresponding to highly radiative regime) to (adiabatic regime), according to our 
equation (O). 
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(iv) We use a differential equation to describe the increase of jet opening angle 9, not 
simply using 9 = a/R. 

Our model also takes many other important factors into account, for example: the equal 
arrival time surface, the distribution of electrons, and the viewing angle. 

We find that during the ultra-relativistic phase, the opening angle 9 increases only 
slightly, but at the Newtonian stage it increases quickly to 9 ~ 0.8 - - 1. The highly 
radiative regime (when e ~ 1) lasts usually for very short period, however a substantial 
fraction of the initial kinetic energy could be dissipated in less than 10 2 - - 10 3 s. As for 
the light curves of optical afterglows, we find no theoretically predicted breaks at the 
critical points where the ejecta transits from 7 > 1/9 to 7 < 1/9, consistent with previous 
numerical studies (Panaitescu & Meszaros 1999; Moderski, Sikora, & Bulik 2000). The 
reason is that at this critical point the lateral expansion just begins to take effect, but its 
action does not prevail over 9q completely, and also the relativistic phase is too short for the 
break to appear. However we have shown that the lateral expansion is even more striking 
during the non-relativistic phase, so we expect an obvious break between the relativistic 
stage and the non-relativistic stage. This has been proved by our calculations made with 
£ e = 0.1. Other calcuations made with £ c = 1 do not show any obvious breaks, but they are 
all characterized by a large timing index. 

We conclude that the most obvious characteristic of beamed eject is the quick decline of 
afterglows at late stages (in fact corresponding to the non-relativistic phase), with a ~ 1.8 
- 2.1, and in some cases, we could even see a sharp break in the light curve between the 
relativistic stage and the non-relativistic stage. We suggest that the quick decay (with a ~ 
1.7 — 2.0) of optical afterglows from some GRBs, such as GRB 970228 (t < 10 d), 980326, 
and 980519, and the breaks in the optical light curves of GRB 990123 and 990510, are 
probably due to beaming effect (Huang, Dai, & Lu 2000a, b, c). 
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But we should also keep in mind that other factors might lead to a large a as well. For 
example, the ISM density may be a decline function of radius (usually n oc R~ 2 ), then an 
isotropic fireball can well explain the quick decay of optical afterglows from GRB 970228, 
980326, and 990519. Additionally, if an isotropic GRB remnant sweeps through a uniform 
ISM (n oc R°) and an uneven ISM (i.e., n oc R~ 2 ) in succession, then an obvious break is 
likely to present in the optical light curve. GRB 990123 and 990510 may also correspond 
to this case. The suggestion by Dai & Lu (1999, 2000) that the quick decay of afterglows 
from GRB 980519 and 990123 is due to an isotropic Newtonian blastwave in a dense ISM 
(with n ~ 10 6 cm -3 ) is also possible. So, the degree of beaming is still very difficult to 
determine now. With the progress in observing technique, when much more GRBs are 
localized, maybe we could also infer some hints from statistical researches. However the 
final solution may come from systematic deep optical surveys, which are expected to find 
many faint decaying optical sources if GRBs are highly collimated. They are afterglows 
from jetted GRBs whose gamma-ray emission deviates the line of sight slightly. 

This work was partly supported by the National Natural Science Foundation of China, 
grants 19773007 and 19825109, the National Climbing Project on Fundamental Researches, 
and the National Project of Fundamental Researches (973 Project). 
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Fig. 1. — Evolution of the Lorentz factor (7). The solid line corresponds to an "ideal" jet 
with "standard" parameters. The dashed line is for "ideal" jet with O — 0.1 and with other 
parameters being equal to those in the "standard" set. For the meaning of "ideal" and 
"standard", please see Section 4 in the main text. 

Fig. 2. — Evolution of the shock radius (R). Parameters and line styles are the same as in 
Fig. 1. 

Fig. 3. — Evolution of the half opening angle (6). Parameters and line styles are the same 
as in Fig. 1. 

Fig. 4. — Schematic illustration of jet evolution on the y — z plane. Z-axis is the symmetry 
axis of the jet and the lateral expansion is approximately at y direction. Parameters and line 
styles are the same as in Fig. 1. 

Fig. 5. — Schematic illustration of the equal arrival time surfaces on the y — z plane. The 
jet is "ideal" and with "standard" parameters. The viewing angle is assumed to be o \> s = 0. 
The dash-dotted lines from right to left correspond to spherical ejecta shells at observing 
time of t — 10 8 s, 10 7 s, 10 6 s, and 10 5 s respectively, enclosed within the jet boundary 
(dotted lines). The solid lines are corresponding equal arrival time surfaces. 

Fig. 6. — Evolution of the radiative efficiency (e) for a "realistic" jet with "standard" 
parameters. 

Fig. 7. — Time dependence of the total kinetic energy (E^) of a "realistic" jet. Parameters 
are the same as in Fig. 6. 

Fig. 8. — R band afterglows. Sr is in units of ergs-s _1 -Hz _1 -cm -2 . The thick solid line 
corresponds to an "ideal" jet with "standard" parameters, and viewing angle 6 ^ s = 0. 
Other lines are drawn with one parameter altered or one condition changed: the thin solid 
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line corresponds to o b S = 0.3; the dashed line corresponds to O — 0.1; the dash-dotted line 
corresponds to £ e = 0.1; and the dotted line is for a "realistic" jet. 
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